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Abstract. I discuss ideal and interacting quantum gases obeying general 
fractional exclusion statistics. For systems with constant density of single-particle 
states, described in the mean field approximation, the entropy depends neither 
on the microscopic exclusion statistics, nor on the interaction. Such systems are 
called thermodynamically equivalent and I show that the microscopic reason for 
this equivalence is a one-to-one correspondence between the excited states of these 
systems. This provides a method, different from the bosonisation technique, to 
transform between systems of different exclusion statistics. In the last section the 
macroscopic aspects of this method are discussed. 

In Appendix A I calculate the fluctuation of the ground state population of 
a condensed Bose gas in grandcanonical ensemble and mean field approximation, 
while in Appendix B I show a situation where although the system exhibits 
fractional exclusion properties on microscopic energy intervals, a rigorous 
calculation of the population of single particle states reveals a condensation 
phenomenon. This also implies a malfunction of the usual and simplified 
calculation technique of the most probable statistical distributions. 
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1. Introduction 

Considerable interest has been shown in the recent years to the study of particle 
systems that exhibit fractional exclusion statistics (FES) - model introduced by 
Haldane in Ref. 0| and which applies, among other, to quasiparticle excitations 
at lowest Landau level in the fractional quantum Hall effect and spinon excitations 
in a spin-i quantum antiferomagnet 0|. The general thermodynamic properties of 
these systems have been deduced mainly by Isakov ^ and Wu [Q. Isakov also 
showed that anyons |^ and one-dimensional (ID) systems of particles described by 
thermodynamic Bethe Ansatz (TBA) 0, ^ have the same thermodynamic behaviour 
as systems exhibiting FES. Since a similar thermodynamic behaviour of two systems 
implies a similarity between the excitation spectra and eventually - depending on the 
ensemble in which the thermodynamics is discussed ~ the same ground state energies 
as functions of the particle numbers, this is irrelevant for comparing the microscopic 
properties of the constituent particles. Sutherland and Iguchi generalized the concept 
of Bethe Ansatz from ID to two and three dimensions (2, 3D), showing that bosons and 
fermions described by this model exhibit also at microscopic level fractional exclusion 
statistics ||]. 

The Bethe Ansatz was first used to calculate hamiltonian eigenvectors of ID spin 
chains with anisotropic spin-spin interaction (see for example Ref. ]lO| and references 
therein) and then applied, starting with Lieb and Liniger Q, to interacting trapped 
particles. The model applies in general to dilute systems, with short-range two-body 
interaction. In the other extreme are dense systems with long range interaction - 
as compared to interparticle distance -, which may be described in the mean field 
approximation (MFA). Nevertheless, MFA is applied also to dilute systems. 

Murthy and Shankar were the first to modify the MFA model for a Fermi system 
[ pl - let me call this new model MFA' - by redistributing the mean field interaction 
among single particle energies, so that a particle of energy e is assumed to interact 



just with particles of energy e' < e (see Appendix B ). Within the MFA' model. 



the interacting Fermi system is equivalent to an ideal Haldane system of statistics 
parameter a which depends on the interaction strength. Independently, Hansson et 
al. jlj] used the MFA' model to define single particle energies in a field theory of 
anyons in the lowest Landau level, while Isakov and Viefers |13| , among other things, 
showed again that this model reproduces Haldane's FES. All these results have been 
obtained for constant DOS in the ideal systems. 

More recently, Bhaduri et al. showed that also a repulsively interacting 2D 
MFA Bose gas has identical thermodynamic parameters as an ideal FES gas with 
a again fixed by the interaction strength. Continuing Ref. jlj], Hansson ct al. 
[p5| discussed the applicability of the MFA to bosons with repulsive delta function 
interaction, in quasi 2D traps, and transformed again the MFA into MFA' to calculate 
the thermodynamic parameters of the corresponding FES gas. 

In what follows I shall say that two systems are thermodynamically equivalent 
if they have the same entropy as a function of temperature, at fixed volume (or 
external potential) and particle number. The amazing thermodynamic similarity in 
2D between interacting Bose or Fermi systems, and ideal FES systems with properly 
chosen statistics parameter a rests on their thermodynamic equivalence. 

To show the special character of nonrelativistic 2D ideal systems, I will start 
by deriving unified expressions for their thermodynamic quantities in terms of 
polylogarithmic functions Km. This bridges Lee's description of bosons and fermions 
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through the intermediate statistics of haldons and emphasizes in a most 
simple way the thermodynamic equivalence of the systems with constant DOS of any 
statistics. The thermodynamic equivalence of 2D Bose and Fermi systems was proved 
by several authors ^ ^ . The fact that the temperature dependence of the 

thermodynamic potentials of Haldane systems with constant DOS is independent of 
statistics was observed before . Here I merely identify the functions involved, which 
makes the writing and manipulations of the thermodynamic quantities much easier 
and compact. In Section 2.2 I give the microscopic explanation of this equivalence, 
which rests not on a microscopic similarity between systems, but on a one-to-one 
correspondence between the excited states of these systems. This also gives us a 
transformation method between systems of different stat istics, which will be extended 
to more general spectra elsewhere In Section 3^, based only on macroscopic 

arguments, I show that all the thermodynamic quantities describing 2D systems in 
the Thomas-Fermi approximation and in any trapping potential, depend only on the 
temperature and the ground state energy density, which is a function of particle 
density. This leads to the definition of classes of thermodynamically equivalent 
systems. Since the proof is based only on macroscopic arguments, it shows explicitly 
that similar thermodynamic behaviour does not imply microscopic similarity. Based 
on Sections ||and||, I can say that the transfo rmation from MFA to MFA' merely sets 
an abstract point of view or counting rule. In Appendix A I calculate the fluctuation 
of the ground state population of a condensed Bose gas in MFA approximation and 
finally, in Appendix B I compare in more detail the results given by MFA and MFA' 
models. I show that in the case of interacting bosons in 2D boxes, depending on the 
choice of a parameter, MFA' may lead to condensation on the lowest energy level, in 
contrast to MFA. 



2. Ideal Haldane systems with constant density of states 
2.1. Thermodynamic equivalence 

In what follows I shall consider Haldane systems with DOS of the form (T(e) = Ce* 
(s > —1), where C is a constant and e is the single particle energy. If the system 
(of nonrelativistic particles) is in a d-dimensional (dD) container with no external 
fields, C is proportional to the hypervolume and s = d/2 — 1. The exclusion statistics, 
characterized by the energy independent parameter a > 0, is manifested between 
particles within the same infinitesimal energy interval Se Q|. The values a — and 
a — 1 correspond to Bose and Fermi statistics, respectively. With these notations, 
the average population of a single particle state is Q 

nie) ^ {wiO + a}-\ (1) 

where w satis fies the equat ion w{(i)"[l + u>(Ce)]"'^^" — (^^ = e'^^^-'/*''^"^, in obvious 
notations (see Appendix B| for a special case) . The grand canonical thermodynamic 



potential, = -PV = keT deCe" log{(l - an)/[l + (1 - a)n]} g, can be put in 
the form 

PV^^ d6 Ce^+'n{e) ^ ^ , (2) 







where U is the internal energy. The total number of particles is iV = Jp°° de Ce*n(e). 
All these functions may be calculated by expressing e in terms of w, but the integrals 
cannot be performed analytically for general s and T. Nevertheless, for s — Q 
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(cr = C), all the thermodynamic quantities can be expressed in terms of elementary 
or poly logarithmic fmictions. I start with 

A^ = fcBTalog(l + 2/o), (3) 

where yo satisfies the equation (l + yo)^~"/yo = = e~'^/^^'^ . Note that yo depends 
on N and T, but not on a. Moreover, = ya{N,T = oo) < yo{N,T) < yo{N,T = 
0) = oo. From Eq. (||) follows 

exp [{^1 - aN/a)/kBT] = 1 - exp [- N / [akYiT)]. (4) 

where I identify the (generalized) Fermi energy as ep = limj^^o M ~ aN / a and observe 
that fjL — e-p IS also independent of a. After some algebra I get the desired expression 
for n and U : 

n = -U = {kBTf(J 



■\o^ {I + yn) + Li2{-yo) 



(5) 

where Li2 is Euler's dilogarithm Using the relation Li2{x) + Li2[—x/{l — a;)] = 
— (1/2) log^ {\ — x), valid for any x <1 p^, one can prove that < for any q > 
and yo > 0, as expected. Combining Eqs. (|^) and (^, it follows: 

^ = -U=^—^— + {kBTfaLi2{-y^), (6) 

In Eqs. (^) and (||) the equivalence between ideal gases obeying any statistics is 
highlighted in the simplest way. Since yo does not depend on statistics, but just on 
N and T, it is obvious that the difference between the thermodynamic potentials of 
particles with different as comes just from an additive constant. All the temperature 
dependence is the same. As an example, one can set a = and a = 1 and use 
the Landen's relations to re-obtain the familiar Bose and Fermi thermodynamic 
potentials [ p^ , |20| . Making use of Eqs. (|) and (|), one can obtain in the usual way 
expressions for the entropy and specific heat, which are both independent of a: 

S = -fc|ra[2iz2(-yo)+log(l + yo)logyo] (7) 

Cv = - - 2klTaLt2hyo). (8) 

Tcr yo 

Since according to Eq. (|^) yo ^ oo when T — > 0, making use of the asymptotic 
behaviour of the dilogarithm, Li2(— yo) ^ + log^(yo)/2] one can show that in 

the limit of low temperatures Cy ~ (7r^/3)fc|Tcr{l — 0[log^(l + yo)/yo]}- In the case 
of particles inside 2D boxes, the leading term is identical to the result for fermions 
obtained by Li et al. |2^. 



2.2. Thermodynamic equivalence from the microscopic point of view: Haldane-Bose 
trans fo rm ation 

Now let me inspect this equivalence from the microscopic point of view. For this I 
consider a Haldane and a Bose system with the same particle number, N, and DOS, 
a = C . In the Haldane system I divide the energy axis into intervals of equal length, 
(5e, and number them, starting from zero, at the lowest interval. Each of these intervals 
contains the same number of single particle states, d = crJe, and a variable number 
of particles, nH,i ( « > 0). For state counting purposes I define a "Bose dimension" 
of the subspace corresponding to the interval i as de.i — d — a(ftH,i — 1), so the 
total size of its Hilbert space is wn.i = {dB,i + "H,i — l)!/['^H,i!(^B,i — 1)!] 0- I 
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Figure 1. Excitation in Bose and Fermi gases with the same equidistant 
single-particle energy levels. The fermions in the g.s., I, II, ...groups are in 
correspondence with the bosons on the ground state, first, second, . . . excited 
states, respectively. In this way I establish a one-to-one correspondence between 
the particle configurations in the Bose and Fermi systems, with the same 
excitation energies. 



denote TVn^i = X)i=o "'H.j and fn.i = ^H,i(* ■ 5e — eF,i), where fn.i shall be called the 
excitation energy of the particle group nn.i, while l-p.i = aN-a^i/a is the (generalized) 
Fermi energy of a similar system, but consisting only of A^H,i particles (I disregarded 
the distribution of particles inside the i'^ energy interval). Obviously, the ground 
state energy of the Haldane system is ?7h,o = YlTLo "H.jCF.i while the total energy 
is C/h = X^i^o "•H.jCH.i — X^i^o^H.i + C%,o- Under canonical conditions, the weight 
of such a distribution is yVH,{rii} = [Hi ^'^P (~'^H,i/fer)i(;H,i]/^H,Af,T, where ^h.w.t 
is the partition function and satisfy the condition X]{n } ^H,{ni} = 1; the sum is 
taken over all distributions {n^}. Now let me put the distributions in the Haldane 
system in correspondence with distributions in the Bose system and calculate the 
new weights. For this, given a distribution {m} in the Haldane system, I divide the 
energy axis of the Bose system in nonequal intervals, 5eB,i = cJB,i/c, and place Ui 
bosons in each of them. As a consequence, the size of the Hilbert space of interval 
i is WB,i — WH,i- Since in a degenerate Haldane gas, in an energy interval 1/ct 
coexist on average 1/a particles, to obtain the correct energy distribution one has 
to overlap a/ a of any consecutive energy intervals above. Then, the excitation energy 

of the group Ui is f^^ X^o ^ ^^^d ~ So ^li^ ~ •^"■j)/"'] = ^H,i- Therefore, the 
two configurations have identical statistical weights. To prove that the two systems 
are equivalent, I have to show that also, to any configuration in the Bose system it 
corresponds one configuration in the Haldane system, with the same excitation energy 
and statistical weight. But this can be done just following the steps outlined above 
in reverse order. In conclusion I showed that there is a one-to-one correspondence 
between the configurations in the Bose and Haldane systems. These configurations 
have the same "excitation energy" and the same statistical weights, so the two systems 
are statistically equivalent even at microscopic level. This analogy finds its simplest 
illustration in the case of ideal Bose and a Fermi gasses with the same spectrum, which 
consist of nondegenerate, equidistant single particle energy levels (ID harmonic trap), 
as shown in Fig. |l| 

From the equivalence proven above and the fact that a macroscopic Bose system 
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with a = C does not experience Bose condensation (see Ref. for interpretations), 
I conclude that at any finite temperature the Fermi sea is not still (excitations occur 
at any depth in the Fermi sea, in macroscopic number). On the other hand, by 
transforming a Bose into a Fermi system, one may describe it (and eventually calculate 
interaction effects) just by considering low energy particle-hole excitations around the 
Fermi surface (bosonisation) , which may considerably simplify the calculations. 



3. Generalization: interacting systems in arbitrary traps 
3.1. Statistical mechanics in boxes 

I investigate here the effects of particle-particle interaction in systems inside d- 
dimensional (dD) containers with no external fields. Next subsection will consider 
arbitrary systems. In the Thomas-Fermi mean field approximation (TF-MFA), the 
particle-particle interaction is replaced by an effective one-body potential, Hi{N) (I 
assume the interaction does not lead to a phase separation) and the total energy of 
the system is 

^a,{k}=2^ek + ^ = iiox + t/^.^.s. + ^ , (9) 

{k} 

where Ck are the single particle energies and Ua^g.s. is the ground-state energy of the 
system without interaction. The quantum numbers that specify the single particle 
states are denoted by k (in the case of spinless particles, k is the momentum). If I 
denote by = ^{k} ^k the energy of the ideal system, then the excitation energy 

(for 2D boxes this is the energy of the equivalent Bose system, as defined in Section 
2.2 ) is i?cx,{k} = -^i°{k} ~ Ua,l.s.- In the microcanonical ensemble N is fixed, so the 
partition function is just a function of the excitation energy Ec^c- -^a,7v(£'Gx)- The 
canonical partition function is usually written as 



Zo..n{I3)^ / dS,e-^^°Z„.^(Sex) (10) 

^^-pNHAN)/2 / d£;(°)e(-''^^"')Z„^A.(i?ex) 



poo 

Jo 



,(/3) 



where I used the shorter notations /3 = (fceT) ^ and Ua,g.s. = Ua/g.s. + ^NHi{N). 
I also denoted by Z'^]^{P) (or Z^'']y(T)) the canonical partition function of the ideal 
Haldane gas of parameter a. Zg°^(/3) corresponds to an equivalent Bose gas of 

excitations. The coefficients in front of Z'^'^]y(/3) and Z'g*'j^(/?) in the two expressions 
for Za^NiP) are redundant in the canonical ensemble, but they should not be ignored 
in grandcanonical calculations. 

Let me now find the population of the single particle energy levels. Since by 
changing the single particle quantum numbers k into k' the total energy of the system 
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changes by ej^. — ek, the relative occupation probabihty of the states k and k' is the 
same as in the case of the ideal Haldane system of parameter a. These occupation 
probabilities in the canonical ensemble follow from z'^]^j{p) which may be calculated 
by the usual saddle point method from the grandcanonical ensemble of the ideal system 
(after dropping the term NHi{N)/2 from the expression of total energy) ~ see Ref. 
[p8| for the special case of Bose-Einstein condensed systems. In this way I obtain the 
average population of the state k, which is given by Eq. dl) with e = and the 
apparent chemical potential, /ia, which in 2D satisfies Eq. (p). Note also that /ia is 
not the real chemical potential, which is defined as 



^ ^ ON dN ON ^ ' 



Taking out the ground state energy from the total energy of the system, as done in the 
expressions ( p^ ) exposes the contribution to the partition function coming just from 
the excitations, which have a bosonic character. Then /ie may be called the "chemical 
potential" of these excitations. Although limj-^o ^o*Ar (7") = 1, for any N (I assumed 
nondegenerate ground state), implies lim^^^o Mb(T') = 0, unlike the chemical potential 
of an ideal Bose system, /zb may take also positive values, as for example in the case of 
ideal ID fermions. In the case of particles in a 2D box, the excitation partition function 



is just the partition function of the equivalent Bose system. In Appendix A I apply 
Eqs. ( pX| ) to calculate the fluctuation of the ground state population in a nonideal 



Bose-Einsten condensed system in grandcanonical ensemble, while in Appendix B I 
discuss in more detail the MFA' model. 



3.2. General trapping potential: thermodynamic considerations 



To apply the procedure introduced in Section 3.1 from a macroscopic point of view, 
I shall first show some general properties of the entropy function of a neutral fluid, 
characterized by the extensive parameters U (internal energy), V (volume), and N 
(particle number): S{U,V, N). Non-essential parts of the proof will be skipped. In 
what follows, the parameters omitted in expressions are the parameters hold fix. 5 is a 
positive, homogenous function of order one, concave downwards |^ [S{XU, XV, XN) = 
XSiU, y, N) and 5(f/i + C/2, Fi + V2,Ni + N2) > 5(C/i, Fi, iVi) + ^(C/a, F2, ^^2)]- The 
inverse function, U{S,V,N) is concave upwards. If I write again the U{S,V,N) = 
Ug.s. (V, N) + Ub (5*, V, N) , then Ug.^. [V, N) = U{S = 0, V, N) is also concave upwards 
(but, in the general case, Ub may not have this property). I take Ug,s.{V,N = 0) = 0. 
Now I fix I assume that A'" and [/(S, Affixed) take values in the intervals [0, 00) 
and [?7g.s.(A^)i 00), respectively (the limits may be restricted further, like for example 
in systems of spins or hard-core particles, but here I work with these intervals). The 
concavity property and the range of U implies that S{U) is monotonic (increasing). 
If I fix t/ = ?7fixcd and t/g.s. = 0, then 5'(C/fixed, N) is also monotonic in A^, otherwise 
is zero at both ends of the allowed interval along the A^ axis. If I introduce the 
function Sb by Sb{Ub,V,N) — S{U,V,N), then Sb and S offer two equivalent 
descriptions of the system. Moreover, if Sb is concave downwards and coincide with 
the entropy of a new system, say B, then we may say that our original system. A, 
is thermodynamically equivalent with B (I denote it by A^ B). In such a case, if 
Ucg.s.iV,N) is any homogenous function of order one, concave upwards, then the 
entropy Sc{Ucg.s.+UB, V, N) = Sb{Ub, V, N) is a legitimate entropy (with all required 
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properties). Moreover, if S'c(t^Cg.s. + Ub,V,N) describes a system C, then A^B^C. 
Therefore for any system described by S{U, V, N) I can define the Bose entropy Sb 
by the procedure outlined above. According to the definition given in Section the 
systems with the same 5b are cahed thermodynamically equivalent and they form an 
equivalence class. 

Since the number of microstates available in a system is dependent on the 
excitation energy Ecx, V, and N [C/g.s. (V, iV) is a redundant quantity], one can 
calculate partition functions having as starting point the Bose entropy of the reservoir, 
iSb- In the grandcanonical ensemble, the probability associated to a microstate 
w{Ecx, V, N) of the system is proportional to exp[fcg^5B(fGx, V, A/")] (I use caUigraphic 
letters for the reservoir quantities). Writing as usual the Taylor expansion of the 
reservoir entropy and identifying Ub = E^x, I obtain the probability of the microstate 
w: 

e-/3(C/B-MBA') 

pH = , (12) 

where Z is the partition function, while /ib is the chemical potential of the bosonic 
excitations, defined by Eq. (pl|). 

Obviously, the 2D Haldane gases described in MFA form an equivalence class. 
Their Bose entropy is the entropy of the equivalent ideal Bose gas. Since /ib is the same 
for all systems in an equivalence class, the difference between the chemical potentials 
of different systems, is due only to the ground state energy (Eq. |ll|). Whether the 
statistics, the interaction, or both, are responsible for the dependence of the ground 
state energy on TV, is of no im portance. Two equivalent gases of different statistics 
(I showed also in Section 3.1 that MFA does not change the microscopic exclusion 
statistics properties), but with the same Ug,s.{N,V), have identical behaviour even 
in general trapping potentials, as long as the Thomas-Fermi approximation holds (so 
that a system in a variable trapping potential can be regarded as a collection of boxes 
with imaginary walls, in contact with each other). 



4. Conclusions 



To summarize, in the beginning of the paper I gave simple expressions, in terms of 
polylogarithmic functions, for the thermodynamic quantities characterizing Haldane 
ideal systems with constant density of single-particle states. These expressions are 
easy to handle and show in a most simple way the thermodynamic equivalence (the 
same entropy function) of such systems. Second, I showed the microscopic reason for 
this equivalence, which is the similarity between the excitation spectra. The model 
used for this purpose is different from the usual bosonisation technique (see Ref. 
for a review) and enables one to transform any Haldane system into a Bose system, and 
vice versa. These transformations: Bose<-^Haldane<-^Fermi, might be very useful if, 
for a given interaction, techniques developed for one kind of systems (like bosonisation 
for Fermi systems) can be transferred to other types of systems (at low temperatures, 
even in Bose systems, the effective thermodynamic contribution comes from a small 
subset of single particle states). 

In the Section ^ I gave a general interpretation of this correspondence technique, 
by subtracting from the internal energy of the system, U, the ground state energy, 
C/g.s. . This led to an equivalent writing of the partition function, in which the internal 
energy is replaced by the excitation energy, U — Ug,s.- In this formulation the chemical 
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potential should be redefined (Eq. |Tl|). Moreover, this redefinition of entropy provides 
a more clear definition of the thermodynamic equivalence, which splits the set of 
physical systems into equivalence classes. Systems belonging to the same class have 
similar excitation spectra (at least up to corrections that vanish in the thermodynamic 
limit). In Appendix A I apply the results of Section ^ to calculate particle fluctuations 



in a non-ideal, condensed Bose system. 



In Appendix B I show some interesting characteristics of what I called the MFA' 
model and discuss similarities and differences with the mean field approximation 
and fractional exclusion statistics. Although on microscopic intervals along the 
pseudoenergy axis, MFA' model reproduces the characteristics of FES gases, as defined 
by Haldane § and Wu [|, in some cases a rigorous calculation of the population 
of single particle states reveals a condensation phenomenon which implies also a 
malfunction of the usual calculation technique. 
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Appendix A. Particle fluctuation in interacting Bose-Einstein gases 



I apply here the technique described in Section || to calculate the fluctuation of 
the ground state occupation number, iVo, of a condensed nonidcal Bose gas, in the 
grandcanonical ensemble. The interaction is assumed to produce a particle dependent 
ground state energy. Let me denote by Z((3, (3^) the grandcanonical partition function. 
Using the notations and arguments from Section ^ I write 



N 



-/3(C/g.,.-MAr) 



> . N 



and I shall omit the subscripts denoting a = 0. Z^'' 
function of the ideal Bose system. As shown in Section |3 



Z)^r' {(3) is the canonical partition 
Ug.s.i^) is a function 



concave upwards. The maximum of the probability distribution over the particle 
number is given by the equation dZ^iP, (3^)/dN = 0, where 



ON 



ZN{l3,l3ti) 



AN 



' /3/i - /3^E 



Z^(/3,/3/i).*Ar(/3,/3/i). 



The solution of the equation will be denoted by iVmax. If there is no iVbound < oo, 
so that limAT^ATi^^^^^^j [dJ/g.s. /dA^] = oo, A^max increases to infinity as [i increases. If 
the system is not pathologic, then hmjv^oo A*b(A^) = and for large enough /x, iVmax 
is the solution of dC/g.s./dA^ — /i. In such a situation we say that the system is 
Bose-Einstein condensed, since for the relevant values of E'cx the configuration space 
does not increase with iV (the ground state is already a particle reservoir). The 
probability distribution Z^r may be expanded around N = iVmax to get Zpf{P,PiJ,) « 
Zw„_(/3,/3Ai) • exp [id^N/dN)N=N^.^AN ~ iVmax)V2] where 



9^jv(/3,/?/i) 
dN 



N=N„ 



diV2 
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For fixed N, the fluctuation of Nq, {S'^Nq)]^, is equal to the fluctuation of the particle 
number on the excited energy levels, {6'^Nex)N (see for example Refs. [|l], |2[ |^ 
for fluctuations in ideal gases). If the probability distribution of Nq in the canonical 
ensemble is P{No,N) = {2Tr{S^ Nq) n)~'^^^ ■ exp[-(No - (Nq) n)^ /2{S^No)n], then the 
grandcanonical probability distribution of A'o is 



2tt dN^ 



2 dN^ 



(A.l) 



Since the system is condensed, the grandcanonical fluctuation of the total particle 
number, (S'^N) is equal to the fluctuation of the canonical average ground state 
occupation number. Using this observation and Eq. (A.l), the total fluctuation of A^o 
in the grandcanonical ensemble follows: 



= / dN 



dNoiNo - {No)n„ 



Therefore the grandcanonic total particle fluctuation of the nonideal gas adds 
quadratically to the canonical fluctuation of the ground state population, to give 
the total mean square fluctuation of the number of particle in the condensate. 



Appendix B. The MFA' model versus mean field approximation and 
fractional exclusion statistics 

If the interaction hamiltonian of Section ^ is Hi{N) = a~^gN, where g is a 
dimensionless constant and a = C (constant), then (mostly) repeating the arguments 
from Refs. ||ll|, [l^ , [ist I define the MFA' model by assimilating the interaction 
energy from Eq. (^ into the single particle (kinetic) energies ek, and I define the 
pseiirfoenergies 

e~(ek) = + V ^/i(e~(ek) - ?(ek'))«(k'), (B.l) 

where rt(k') is the population of the state k' and h{x) — 1, ho, or 0, depending 
weather x > 0, a; = 0, or a; < 0, respectively; in what follows I assume g > —a 
and ho G [0,1], where a is the statistics parameter of gas without interaction.^ As 
Wu I divide the energy and pseudoenergy axes into infinitesimal (microscopic) 
intervals. If exist such microscopic divisions for which none of the intervals is 
completely blocked, then the energy interval [et, ek + de) is transformed into the 
pseudoenergy interval [ek, ek + de), where de = de + a~'^gp{e-k}de = de + a^^ gp{e)de, 
and (Tn(ek)de = p(ek)de = p(e)de represent the average number of particles in the 
intervals under consideration. Moreover, the Bose dimensions of the two intervals de 
and de are: 

^B.de = crde - Q;p(ek)de = crde - (a + g)p{e)de = dB.de, (B.2) 

where a — gp{e) is the effective density of states along the pseudoenergy axis. The 
Bose density d-Q.de is identical to the Bose density of an ideal Haldane gas of parameter 

I ^In the counting rule of Hansson et al. , two particles could not have the same pseudoenergy 
even if they occupy the same single particle state. Since I workwith a (quasi)continuous spectrum, 
I adopt the more general point of view of Isakov and Viefers [13[ , but relaxing also the constraint 
ho = \. 
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a + g (see also Ref. |Q for related discussions). This implies, in accordance with Refs. 
[l3|, that the MFA' model leads to the same microscopic particle distribution /5(e) 



as an ideal Haldane gas of statistics parameter a + g, call it na+g{e) (|l|) (I included the 
subscript a + g in the notation to specify the statistics parameter of an ideal Haldane 
gas). 

Appendix B.l. 2D Bose gas in MFA' approximation: an unexpected condensation 

To continue exercising with the MFA' approximation, let me now consider a 2D 
interacting Bose gas. I assume that the particle-particle interaction is repulsive, 
so Hi{N) = a~^gN, with g > 0. Under canonical condition and in the MFA 
approximation, the gas is thermodynamically equivalent with the ideal Bose gas and 
its entropy and specific heat are given by the universal expressions (|^) and (||). On 
the other hand, in the MFA' model the energies of the single particle states are 



"shifted" according to Eq. (B.l), but the occupation of the single particle states 
is done in accordance with the "original" Bose statistics (see Eq. |B.2| ). If I number 
the single particle (kinetic) energies and pseudoenergies as eg = < ei < . . ., and 
Co < £i < ■ • respectively (I assume level is nondegenerate) , then the MFA' 
grandcanonical partition function is 

{m} i>l 

no 

where zix'^ = X]{n } ni>i > while Ui represents the population of the 

single particle state i. The exponent —Pa^^ghoUQ + f3^nQ has a maximum for 
no = no,max = ■ (25/10)"^ ■ If < 0, no, max < 0, and I expect no macroscopic 



population of any single particle state. In such a case, according to Eq. ( |B.2D , 
p{e) — ang{e) and I should recover the previous results for the ideal Haldane system 
of parameter g (see Q). This will be shown to be true and for the Haldane gas of 
parameter g > (with the chemical potential /Lin) under canonical conditions, exists a 
strictly positive temperature Tinv (defined by Eq. ^), so that Hg{T < Tinv) > 0. Since 
for T > Tinv, M = Mh, I conclude that /i(T > Tjnv) ^ as T ^ Tinv Turning 
back to the MFA' model in grandcanonical ensemble, I investigate the situation 
/X > 0. li fj, > 0, the distribution e"'^'^ ghgiig+Pfma j^^^g ^ very sharp maximum 
centred at no — no, max > and with the second moment = (/3^)~^no,max- The 
relative fluctuation of no is \/<5^ no/no, max = (/3m)~^''^''^o max' which is zero in the 
thermodynamic limit. Most often, such a sharp distribution is believed to fix the 
average value of the variable at the distribution maximum, but this is an example 
where this general procedure would lead to wrong results. 

For a rigorous calculation, I denote by Zg{(},(3^) the grandcanonical partition 
function of the Haldane system of parameter g and I assume for the beginning that 
ni>i are all microscopic populations (in the end this will turn out to be true). Then 

Z{P,p^i) ^ ^i„,(/3,/3Ai) =^e-^'^"'5'"'"o+^'^"° •Z,[/3,/3(a*-?i)] 

no no 
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The most probable value of uq in the summation (B.4) is given by the equation 
dZng/duf) = 0. To find this value I define the function 



h,{no, = Z-^^ = ~2hono'-^ + /?M - (A^cx)^, (B.5) 



-i^Zno „, fig , „ n^T \l^9 

— = -2/iono — ' '^^ ^ 

ano a 

where (A^ox) = d\og{Zg)/d{(3fi') is the average number of particles on the excited 
energy levels. The most probable value of no is given by one of the zeros of /. For 
simplicity I shall use the notations ^ = f3nQg/a, C, = and ^' = ji — a^^gno- 
Plugging Eq. (|) into I can rewrite / in two equivalent ways: 



ho (C, C) = log [jf^J + (1 - 2/io)e (B.6) 
= (1 - 2ho)C + 2ho \og{y'o) - [2ho + g{l - 2ho)] log(l + y^,), (B.7) 
where (1 + v'qY~^ ly'o = e~'^^~^\ Expression (B^) imphes that fhoi^X) < for any 



ho and fhoi^X) < for any ^, if ho > 1/2. Therefore the probability distribution of 
no has a local maximum at riQ — 0. If Hq > 1/2, then rip = is the only maximum of 
fho- Since Z„„ ^ Zg(/3, /3^)-exp (^g-^ksTa J^^''"^ //io(C',C)d^') and ksTa is assumed 

to be a macroscopic quantity, it follows that Zng has an infinitely sharp maximum at 
no = which implies that the ground state is microscopically populated in spite of 
the very sharp maximum of e~^'^ ghoug+iStma ^ centred at no. max- 

To observe the behaviour of fhg for large f we have to look at the expression 



(B.7). For ho = 0, the function gets the simple form /o(CjC) = C ~ 5log(l + y'o)- 
Since yo(C ~ decreases monotonically from y'oiC) to ?/o(~oo) = 0, as ^ increases 
from to oo, then /o(C,C) is also a monotonic function of ^, with /o(0, C) < and 
/o(oo, ^) = ^ > 0. The only zero of fo is at j/q = y'^ = e^/^ — 1 and corresponds to 



a minimum of probability. As one can observe directly from Eq. (B.3), the maximum 
probability is cxi and corresponds to np = cxd. Yet, as I mentioned above, the system 
has a metastable state for microscopic no, which corresponds to the equivalent (ideal) 
Haldane distribution. 

The function f^o {£,, C) is a continuous function in all the parameters and variables. 
From Eq. ( p3.6[) fo llows that for fixed ^ and / is a decreasing function of ho 
and from Eq^ ( |B.7| ) we observe that for '3> C: fhg C) ~ C ~ 2ho£,. For small 
enough ho, fhoitO = has a solution at ^max ~ C/{2ho) and exp(C - Cmax) = 
exp {— [(2/io)~^ ~ C} ^ exp(C/g) — 1. In such a case, and since fho>o{£, ^ oo, C) ^ 
—00 and also fho (0, C) < 01 can conclude that / has two zeros. The first corresponds to 
a local minimum of probability for no , while the second to a local maximum. Whether 
no = or no = no^max = Cmaxf/ (/3<?) ~ (^f^\2gho)~^ has higher probability, it depends 
on the specific values of ho and C- Nevertheless, continuity of / and monotonicity in 
ho implies that exists a critical value of /iq, between and 1/2 and which increases 
monotonically with C, I call it ft,o,cr(/3M)i so that for ho < /io,cr(/3M) the ground state 
is macroscopically populated (the probability for no — no, max is highest), while for 
ho > ^o,cr the ground state has microscopic occupation number and is well described 
by the Haldane ideal gas of parameter g. If the ground state is macroscopically 
populated, I say that the 2D MFA ' Bose system is condensed. 

Now I prove that none of the excited energy levels are macroscopically occupied. 
Obviously, I start with ei. This level can be macroscopically occupied if and only if 
A*'("-o,max) > and ho < ho.aiP^i')- The condition /i' > implies ho < ho.aiPfj) and 
'V^. But for no = ^ = C, and fhoiCO = log[2/o(0)] - log[l + 
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2/o(0)] + (1 — 2/io)C - where yo{0) > 1 is a fixed value which satisfy the equation 
2/o(0) • (1 + yo(0))^^^ = 1- Therefore fhoiC: C) is a function Hnear in which starts at 
fhoiO, 0) = log[yo(0)] -log[l + ?/o(0)] < and, since ho < ho,cAPfA < 1/2, increases to 
infinity, as C increases. The continuity of fhg {£,, C) completes then the proof that always 
C < Cmax, which implies that no, max > <7~^cr/i. This proves the fact that ei and as a 
consequence any state of pseudoenergy ei>i is microscopically populated. Also, the 
monotonic increase of ft.o,cr with C implies that, if the ground state is microscopically 
populated, so are all the other states. 

To finish the exercise I will show that for any < ho < 1/2, exists a temperature 
Tc < Tinv at which the system condenses on the ground state. If the system is not 
condensed, it is described as a Haldane gas of parameter g. In this case the relative 
fluctuation of the total particle number N vanishes in the thermodynamic limit. On 
the other hand, for a condensed system, 

2no = • exp I — — J fh„ , C)d^ I 



2^no.max • exp I —("0 - no^max) 



2a ok, 



(B., 



where I used the fact that //io(^max, C) — 0- From (|B.8| ) it is easy to observe that the 

— 1/2 

relative fluctuation is proportional to rip J^^^, which vanishes in the thermodynamic 
limit. The relative total particle fluctuation, which is obtained by adding together the 
contributions from the ground state, from the excited states (described as a Haldane 
gas), and the correlations between them, vanishes also. It is therefore safe to use 
the grandcanonical average values in canonical calculations even for this unusual toy 
system (fine-tuning due to the change of ensemble are relevant only for the finite 



size effects). At T — Tinv, fJ- = and so is C- In this case, from Eq. (B.7) we 
see that fho{(,,0) < for any ^ and ho, therefore no is still microscopic, so the 
condensation temperature is lower. If the system does not condense, /i increases 
as the temperature is lowered, which implies an even faster increase of C- If I assume 
that for a chosen ho < 1/2, — (the system does not condense), then I can choose 
T so that C takes any value up to infinity. If C — ^ ^ 1, then y'^ ~ eiC^O/g and 
log[(l + y'o)/y'o\ ~ e"(^~^)/s < 1. Now, for any ho < 1/2 and interval [0,^o], I can 
choose C » 1 + Co, so that fh,X^o, C) > and log[(l + y'o)/y'o] < (1 - 2ho)^^/2, 
which implies that -Zno.max > ^no(Co) > = 2,g{l3,f3^). In this case the system 
is condensed, with ng^max particles on the ground state, so the initial assumption 
(Tc = 0) was wrong. 

In conclusion, for any ho < 1/2, exists a temperature Tc € (0,Tinv), at which 
the condensation occurs and bellow which the similarity between the MFA' Bose 
system and the "usual" ideal Haldane system Q of parameter g is lost. Moreover, the 
onset of the condensation removes also any thermodynamic equivalence between gases 
described in MFA and MFA' models. The exercise presented here is also interesting 
for the fact that it showed in a concrete example how an infinitely sharp (in the 
thermodynamic sense) probability distribution of particle on the ground state may be 
overwhelmed by the probability distribution of particles on the excited energy levels. 
Vice versa, it also shows that the equilibrium distribution of the Haldane gas, as it 
was deduced by Isakov H, |j and Wu [Q, and which corresponds to microscopic no. 



may not be the equilibrium distribution, in spite of the FES manifested in Eq. (B.2) 
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